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For states in infinite dimensional Hilbert spaces entanglement quantities like the entanglement of 
distillation can become infinite. This leads naturally to the question, whether one system in such 
an infinitely entangled state can serve as a resource for tasks like the teleportation of arbitrarily 
many qubits. We show that appropriate states cannot be obtained by density operators in an 
infinite dimensional Hilbert space. However, using techniques for the description of infinitely many 
degrees of freedom from field theory and statistical mechanics, such states can nevertheless be 
constructed rigorously. We explore two related possibilities, namely an extended notion of algebras 
of observables, and the use of singular states on the algebra of bounded operators. As applications 
we construct the essentially unique infinite analogue of maximally entangled states, and the singular 
state used heuristically in the fundamental paper of Einstein, Rosen and Podolsky. 

Dedicated to the memory of Rob Clifton 



I. INTRODUCTION 

Many of the concepts of entanglement theory were orig- 
inally developed for quantum systems described in finite 
dimensional Hilbert spaces. This restriction is often jus- 
tified, since we are usually only trying to coherently ma- 
nipulate a small part of the system. On the other hand, 
a full description of almost any system, beginning with 
a single elementary particle, requires an infinite dimen- 
sional Hilbert space. Hence if one wants to discuss de- 
coherence mechanisms arising from the coupling of the 
"qubit part" of the system with the remaining degrees of 
freedom, it is necessary to widen the context of entangle- 
ment theory to infinite dimensions. This is not difficult, 
since many of the basic notions, e.g. the definitions of 
entanglement measures, like the reduced von Neumann 
entropy or entanglement of formation, carry over almost 
unchanged, merely with finite sums replaced by infinite 
series. More serious are some technical problems arising 
from the fact that such entanglement measures can now 
become infinite, and are no longer continuous functions 
of the state. Luckily, as shown in recent work of Eisert 
et. al. 0], these problems can be tamed to a high degree, 
if one imposes some natural energy constraints on the 
systems. 

In the present paper we look at some not-so-tame 
states, which should be considered as idealized descrip- 
tions of situations in which very much entanglement is 
available. For example, in the study of "entanglement 
assisted capacity" |2| one assumes that the communicat- 
ing partners have an unlimited supply of shared maxi- 
mally entangled singlets. In quantum information prob- 
lems involving canonical variables it is easily seen that 
perfect operations can only be expected in the limit of 
an "infinitely squeezed" two mode gaussian state as en- 
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tanglement resource (see also Section VI). But infinite 
entanglement is not only a desirable resource, it is also 
a natural property of some physical systems, such as the 
vacuum in quantum field theory (see ||, |j and Section |y| 
below). Our aim is to show that one can analyze these 
situations by writing down bona fide states on suitably 
constructed systems. 

It turns out that in order to do this we need to go 
one step beyond the standard Hilbert space formalism 
of quantum mechanics (even with infinite dimensional 
Hilbert spaces). This is completely analogous to other 
systems with infinitely many degrees of freedom, aris- 
ing in quantum statistical mechanics (in the thermody- 
namic limit) and quantum field theory. For example, 
consider a translationally invariant, finite density equi- 
librium state of a gas. For such a state the probability 
for having only finitely many particles must obviously 
be zero. In contrast, a state given by density operators 
on Fock space implies a probability distribution of parti- 
cle numbers, which are by definition finite numbers in N. 
The expectation value of this distribution may be infinite. 
This means that in the limit of taking many (m — > oo) 
equally prepared systems, the total number n of particles 
diverges faster than m, so that n/m — > oo. However, the 
particle number of each single system remains finite. 

Our first main result establishes a similar limitation 
for entanglement theory. Consider a state with infinite 
distillible entanglement. This means that when we take a 
large number m of copies of such system pairs, and try to 
obtain from these a large number n of nearly maximally 
entangled singlet pairs by local operations and classical 
communication, we can achieve a rate n/m-xxiasm-* 
oo. In this asymptotic sense each pair contains infinitely 
many ebits of entanglement. But what does this mean for 
a single pair? Should we not get arbitrarily many singlets 
even out of this? According to Theorem [l] the answer is 
no, as long as we stay in standard quantum mechanics 
with bipartite states given by density operators on tensor 
product Hilbert spaces (of finite or infinite dimension). 

However, if we follow the lead of statistical mechan- 
ics, and employ the methods for describing states with 
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actually infinite particle number, we find a very natu- 
ral framework to overcome this limitation: Here we can 
have a direct mathematical representation of the intu- 
itive idea of "having infinitely many singlets" . But we 
need to reconsider either the notion of states, allowing 
"singular" probability assignments on the space of quan- 
tum observables, which canno t be written in terms of 
density operators (cf. Section [V B ) or else allow more 
gene ral o bservable algebras for Alice and Bob (cf. Sec- 
tion IV C). These two approaches are closely related. We 



show in Section how this extended framework leads 
to an essentially unique description of maximally entan- 
gled states of systems with infinitely many degrees of 
freedom. In the final Section VI we discuss the "origi- 
nal EPR-state", i.e., a mathematically rigorous version 
of the singular state employed by Einstein, Rosen and 
Podolsky in their fundamental paper H. This Section 
builds on the work of Rob Clifton || ]?], |sj , who sadly 
died while this paper was in preparation. We dedicate it 
to his memory. 



II. DENSITY OPERATORS ON INFINITE 
DIMENSIONAL HILBERT SPACE 

We will start our discussion with a short look at en- 
tanglement properties of density operators on an infinite 
dimensional but separable 1 Hilbert space H^Ti.. 

Most of the definitions of entanglement quantities 
carry over from the finite dimensional setting without 
essential change. Since we want to see how these quan- 
tities may diverge, let us look mainly at the smallest, 
the distillible entanglement. It is defined as the largest 
rate n/m of n nearly perfect singlets, which can be ex- 
tracted from m pairs prepared in the given state by pro- 
tocols involving local quantum operations and classical 
communication. The class of quantum operations is de- 
fined as in the finite dimensional case, by completely pos- 
itive, trace preserving maps (in the Schrodinger picture). 
Effectively, distillation protocols for infinite dimensional 
systems can be built up by first projecting to a suitable 
finite dimensional subspace, and subsequently applying 
finite dimensional procedures to the result. 

With this in mind we can easily construct pure states 
with infinite distillible entanglement. Let us consider vec- 
tors in Schmidt form, i.e., 



with orthonormal bases e' 



> 0, and J2 n 



(1) 



and positive numbers 



1. The density operator of the 



1 Another extension of this framework, namely to Hilbert spaces 
of uncountable dimension (i.e., unseparable ones in the topolog- 
ical sense) is not really interesting with regard to entanglement 
theory, since any density operator has separable support, i.e., it 
is zero on all but countably many dimensions. 



restriction of this state to Alices's subsystem has eigen- 
values c 2 , and von Neumann entropy —^2 n c^ l \og 2 (c^ l ), 
which we can take to be infinite (c„ = l/(Z(n + 
2)log 2 (rt + 2) 2 ) will do). We can distill this by using 
more and more of the dimensions as labelled by the bases 
e' n ,e'n, and applying the known finite dimensional distil- 
lation procedures to this to get out arbitrary amount of 
entanglement per pair. 

Once this is done, it is also easy to construct mixed 
states with large entanglement in the neighborhood of 
any state p, mixed or pure, separable or not. We only 
have to remember that every state is essentially (i.e., up 
to errors of small probability) supported on a finite di- 
mensional subspace. Therefore we can consider the mix- 
ture p e = (1 — e)p+eo~ with a small fraction of an infinitely 
entangled pure state er, which is supported on those parts 
of Hilbert space, where p is nearly zero. Therefore dis- 
tillation based on the support of a will work for p e and 
produce arbitrarily large entanglement per p e pair, in 
spite of the constant reduction factor e. 

For the details of such arguments we refer to |t|. [|. 
The argument as given here does not quite show that 
states of infinite distillible entanglement are norm dense, 
but it certainly establishes the discontinuity of the func- 
tion "distillible entanglement" with respect to the trace 
norm topology. This might appear to show that the ap- 
proach to distillible entanglement based on finite dimen- 
sional systems is fundamentally flawed: If only finitely 
many dimensions out of the infinitely many providing a 
full description of the particle/system are used, might not 
the entanglement be misrepresented completely? Here it 
helps that states living on a far out subspace in Hilbert 
space usually also have large or infinite energy. For typi- 
cal confined systems, the subspaces with bounded energy 
are finite dimensional, so if we assume a realistic a pri- 
ori bound on the energy expectation of the states on the 
consideration, continuity can be restored [EJ. 



III. INFINITE ONE-COPY ENTANGLEMENT 

If the entanglement of formation of a state is infinite: 
how much of that entanglement can we get out? Since 
for pure states the distillible entanglement is the same 
as the entanglement of formation we know that given 
sufficiently many copies of the state, we can use a distil- 
lation process producing in the long run infinitely many 
nearly pure singlets per original entangled pair. But if 
the entanglement is infinite, might it not be possible to 
use only one copy of the state in the first place? In other 
words, are there states, which can be used as a one time 
resource, to teleport an arbitrary number of qubits? 

We will now give a definition of such states. The ex- 
traction of entanglement will be described by a sequence 
of operations resulting in a pair of d-level systems with 
finite d. The extraction is successful, if this pair is in a 
nearly maximally entangled state, when one starts from 
the given input state. The overall operation is then given 
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mathematically by a completely positive, trace preserv- 
ing map Ed- Of course, we must make sure that the ex- 
traction process does not generate entanglement. There 
are different ways of expressing this mathematically. For 
example, we could allow £d to be composed of an arbi- 
trary number of rounds of local quantum operations and 
classical communication ("LOCC operations"). We will 
also consider a much weaker, and much more easily veri- 
fied condition, namely that £d takes pure product states 
into states with positive partial transpose ("PPtPPT 
operations" for "pure product to positive partial trans- 
pose"). Of course, every LOCC channel is a PPtPPT 
channel. 

The success is measured by the fidelity (overlap) of the 
output state £d{p) with a fixed maximally entangled state 
on C d ® C d . By pd we denote the projection onto this 
maximally entangled reference vector. Then a density 
operator p is said to have "infinite one-copy entangle- 
ment" , if for any e > and any d £ N there is a PPtPPT 
channel £d such that 

tr{£d( P )pd) > 1 - e . (2) 

Then we have the following Theorem, whose proof uses a 
distillation estimate of Rains |hJ developed for the finite 
dimensional context. 

Theorem 1 For any sequence of PPtPPT channels £d, 
d £ N, and for any fixed density operator p we have 

km tr{£ d {p)pd) = . (3) 

d — >oo 

In particular, no density operator with infinite one-copy 
entanglement exists. 

Proof. Consider the operators Ad defined by 

tx{pA d ) = tx{£ d {p)p d ). (4) 

In order to verify that A d exists, observe that £d, as a pos- 
itive operator is automatically norm continuous. Hence 
the right hand side is a norm continuous linear functional 
on density matrices p. Since the set of bounded opera- 
tors is the dual Banach space of the set of trace class 
operators |ll], Theorem VI. 26] such functionals are in- 
deed of the form (|J). We now have to show that, for 
every p, we have lim^ tx(pAd) = 0, i.e., that Ad — > in 
the weak*-topology of this dual Banach space. 

Obviously, < A d < 1, and by the Banach- Alaoglu 
Theorem JO], Theorem IV. 21], this set is compact in 
the topology for which we want to show convergence. 
Hence the sequence has accumulation points and we only 
have to show that all accumulation points are zero. Let 
Aac denote such a point. Then it suffices to show that 
tx(aA oc ) = for all pure product states a. Indeed, since 
Aoo > 0, this condition forces Aoo^xg) if> = for all pairs 
of vectors 0, ip, and hence A m = 0, because such vectors 
span the tensor product Hilbert space. 



On the other hand, our locality condition is strong 
enough to allow us to compute the limit directly for pure 
product states a. We claim that 

tx(aA d ) = tr(£d(cr)pd) 

= tx(£ d (af 2 p^)<\\pf\\=l/d (5) 

Here we denote by X T2 the partial transposition with 
respect to the second tensor factor, of an operator X on 
the finite dimensional space C d ® C d , and use that this 
operation is unitary with respect to the Hilbert-Schmidt 
scalar product (X,Y)hs = tx(X*Y). By assumption, 
£d{o) T2 > 0, and since partial transposition preserves 
the trace, £d{c) T2 is even a density operator. Hence the 
expectation value of p T d 2 in this state is bounded by the 
norm of this operator. But it is easily verified that p T d 2 
is just (1/d) times the unitary operator exchanging the 
two tensor factors. Hence its norm is {1/d). Taking the 
limit of this estimate along a sub-net of Ad converging to 
Aoo, we find tr(crA oc ) =0. □ 



IV. SINGULAR STATES AND INFINITELY 
MANY DEGREES OF FREEDOM 

In this section we will show how to construct a sys- 
tem of infinitely many singlets. It is clear from Theo- 
rem [j] that not all of the well-known features of the finite 
situation will carry over. Nevertheless, we will stay as 
closely as possible to the standard constructions trying 
to pretend that oo is finite, and work out the necessary 
modifications as we go along. 

A. Von Neumann's incomplete infinite tensor 
product of Hilbert spaces 

The first difficulty we encounter is the construction of 
Hilbert spaces for Alice's and Bob's subsystem, respec- 
tively, which should be the infinite tensor power (C 2 )®°° 
of the one qubit space C 2 . Let us recall the definition of a 
tensor product: it is a Hilbert space generated by linear 
combination and norm limits from basic vectors written 
as $ = 0jj where (j>j is a vector in the jth ten- 

sor factor. All we need to know to construct the tensor 
product as the completion of formal linear combinations 
of such vectors are their scalar products, which are, by 
definition, 

/ oo oo \ oo 

\j'=l J=l / J=l 

The problem lies in this infinite product, which clearly 
need not converge for arbitrary choice of vectors <f>j , ipj . 
A well-known way out of this dilemma, known as von 
Neumann's incomplete tensor product ]l2| is to restrict 
the possible sequences of vectors <f>\ , 4>2 , ■ • ■ in the basic 
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product vectors: for each tensor factor, one picks a ref- 
erence unit vector \j > and only sequences are allowed for 
which cfij — Xj holds for all but a finite number of indices. 
Evidently, if this property holds for both the <j>j and the 
tpj the product in (Q) contains only a finite number of 
factors 7^ 1, and converges. By taking norm limits of 
such vectors we see that also product vectors for which 
J27Li 1 1 "A? — XjII < 00 are included in the infinite product 
Hilbert space. However, the choice of reference vectors 
Xj necessarily breaks the full unitary symmetry of the 
factors, as far as asymptotic properties for j — ► 00 are 
concerned. For the case at hand, i.e., qubit systems, let 
us choose, for definitcness, the "spin up" vector as Xj f° r 
every j, and denote the resulting space by Hoo- 

An important observation about this construction is 
that all observablcs of finite tensor product subsystems 
act as operators on this infinite tensor product space. In 
fact, any operator <S)^Li Az makes sense on the incom- 
plete tensor product, as long as Aj = 11 for all but finitely 
many indices. The algebra of such operators is known as 
the algebra of local observables. It has the structure of 
a *-algebra, and its closure in operator norm is called 
quasi-local algebra fl3| . 

Let us take the space Hoo as Alice's and Bob's Hilbert 
space. Then each of them holds infinitely many qubits, 
and we can discuss the entanglement contained in a den- 
sity operator on Hoo <8> Hoo ■ Clearly, there is no general 
upper bound to this entanglement, since we can take a 
maximally entangled state on the first M < 00 factors, 
complemented by infinitely many spin-up product states 
on the remaining qubit pairs. But for any fixed den- 
sity operator the entanglement is limited: for measure- 
ments on qubit pairs with sufficiently large j we always 
get nearly the same expectations as for two uncorrelated 
spin-up qubits (or whatever the reference states Xj dic- 
tate). This is just another instance of Theorem [l| there 
is no density operator describing infinitely many singlets. 



B. Singular states 

However, can we not take the limit of states with grow- 
ing entanglement? To be specific, let $m denote the 
vector which is a product of singlet states for the first 
M qubit pairs, and a spin-up product for the remaining 
ones. These vectors do not converge in Hoo <8 Hoo , but 
that need not concern us, if we are only interested in ex- 
pectation values: for all local observables A (observables 
depending on only finitely many qubits), the limit 

v(A) = lim($M,A$A/) (?) 

exists. Thereby we get an expectation value functional 
for all quasi-local observables, and by the Hahn-Banach 
Theorem (see e.g. [Ell Theorem III.6]), we can extend 
this expectation value functional to all bounded opera- 
tors on Hoo ® Hoo ■ The extended functional u> has all 
the properties required by the statistical interpretation 



of quantum mechanics: linearity in A, oj(A) > for pos- 
itive A, and uj(t) = 1. In the terminology of the theory 
of operator algebras, it is a state on the algebra of all 
bounded operators. By construction, u> describes maxi- 
mal entanglement for any finite collection of qubit pairs, 
so it is truly a state of infinitely many singlets. 

How does this match with Theorem |l|? The crucial 
point is that that Theorem only speaks of states given 
by the trace with a density operator, i.e., of functionals 
of the form u> p (A) = ti(pA). Such states are called "nor- 
mal". But there is no density operator for lo: this is a 
singular state on the algebra of bounded operators. 

Singular states are not that unusual in quantum me- 
chanics, although they can only be "constructed" by an 
invocation the Axiom of Choice, usually through the 
Hahn-Banach Theorem 2 . For example, we can think of 
a non-relativistic particle localized at a sharp point, as 
witnessed by the expectations of all continuous functions 
of position. Extending from this algebra to all bounded 
operators, we get a singular state with sharp position 3 , 
but "infinite momentum", i.e., the probability assigned 
to finding the momentum in any given finite interval is 
zero fHI . This shows that the probability measure on the 
momentum space induced by such a state is only finitely 
additive, but not er-additive. This is typical for singular 
states. 

More practical situations involving singular states arise 
in all systems with infinitely many degrees of freedom, as 
in quantum field theory and in statistical mechanics in 
the thermodynamic limit. For example, the equilibrium 
state of a free Bose gas in infinite space at finite density 
and temperature is singular with respect to Fock space 
because the probability for finding only a finite number 
of particles in such a state is zero. In all these cases, 
one is primarily interested in the expectations of certain 
meaningful observables (e.g., local observables), and the 
wilder aspects of singular states are connected only to the 
extension of the state to all bounded operators. There- 
fore it is a good strategy to focus on the state as an 
expectation functional only on the "good" observables. 



C. Local observable algebras 

If we want to represent a situation with infinitely many 
singlets, an obvious approach is to take again von Neu- 
mann's incomplete tensor product, but this time the in- 
finite tensor product of pairs rather than single qubits, 
with the singlet vector chosen as the reference vector Xj 



2 Other constructions based on the Axiom of Choice are the ap- 
plication of invariant means, e.g., when averaging expectation 
values over all translations, or algebraic constructions using max- 
imal ideals. For an application in von Neumann style measure- 
ment theory of continuous spectra, see |l4| 

3 This is not related to improper eigenkets of position, which do 
not yield normalized states 
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for every pair. We denote this space by Hoooo, and by 
G TLoaoa the infinite tensor product of singlet vectors. 
Clearly, this is a normal state (with density operator 
|0)(0|), and we seem to have gotten around Theorem [j] 
after all. 

However, the problem is now to identify the Hilbert 
spaces of Alice and Bob as tensor factors of 7ioooo- To 
be sure, the observables measurable by Alice and Bob, 
respectively, are easily identified. For example, the a x - 
Pauli matrix for Alice's 137 th particle is a well defined 
operator on Hocoa- Alice's observable algebra A is gen- 
erated by the collection of all Alice observables for each 
pair. Bob's observable algebra B is generated simi- 
larly, and together they generate the local algebra of 
the pair system. Moreover, the two observable alge- 
bras commute elementwise. This is just what we expect 
from the usual setup, when the total Hilbert space is 
H = Ha ®Hb, and Alice's and Bob's observable alge- 
bras are A = B(Ha) ® Ifl and B = 1 A <8> B(H B )- 

However, the A and B constructed above are definitely 
not of this form, so 7ioooo has no corresponding decom- 
position as Ha ® Hb ■ The most direct way of seeing this 
is to note that Hoooo contains no product vectors describ- 
ing an uncorrelated preparation of the two subsystems. If 
we move to qubit pairs with sufficiently high index, then 
by construction of the incomplete tensor product, every 
vector in 7ioooo will be close to the singlet vector, and in 
particular, will vi olat e Bell's inequality nearly maximally 
(see also Section VB). 

Hence we arrive at the following generalized notion of 
bipartite states, generalizing the finite dimensional one: 
Alice's and Bob's subsystems are identified by their re- 
spective observable algebras A and B. We postpone the 
discussion of the precise technical properties of these al- 
gebras. What is important is, on the one hand, that these 
algebras are part of a larger system, so they are both sub- 
algebras of a larger algebra, typically the algebra B(H) 
of bounded operators on some Hilbert space. This al- 
lows us to consider products and correlations between 
the two algebras. On the other hand, each measurement 
Alice chooses must be compatible with each one chosen 
by Bob. This requires that A and B commute elemen- 
twise. A bipartite state is then simply a state on the 
algebra containing both A and B. 

We can then describe the two ways out of the NoGo- 
Theorem: on the one hand we can allow more general 
states than density matrices, but on the other hand we 
can also consider more general observable algebras. In 
the examples we will discuss, the algebra containing A 
and B will in fact be of the form B(H), and the states 
will be given by density matrices onH. So both strategies 
can be successful by themselves. 



D. Some basic facts about operator algebras 

The possibility of going either to singular states or to 
extended observable algebras is typical of the duality of 



states and observables in quantum mechanics. There are 
many contexts, where it is useful to extend either the 
set of states or the set of observables by idealized ele- 
ments, usually obtained by some limit. However, these 
two idealizations may not be compatible Jl5| . There are 
two types of operator algebras which differ precisely in 
the strength of the limit procedures under which they are 
closed . 

On the one hand there are C*-algebras, which are iso- 
morphic to norm and adjoint closed algebras of operators 
on a Hilbert space. Norm limits are quite restrictive, so 
some operations are not possible in this framework. In 
particular, the spectral projections of an hermitian ele- 
ment of the algebra often do not lie again in the algebra 
(although all continuous functions will). Therefore, it is 
often useful to extend the algebra by all elements ob- 
tained as weak limits (meaning that all matrix elements 
converge). In such von Neumann algebras the spectral 
theorem holds. Moreover, the limit of an increasing but 
bounded sequence of elements always converges in the 
algebra. For these algebras the distinction between nor- 
mal and singular states becomes relevant. The normal 
states are simply those for which such increasing limits 
converge, and at the same time those which can be rep- 
resented by a density operator in the ambient Hilbert 
space. 

A basic operation for von Neumann algebras is the for- 
mation of the commutant: for any set M. C B(H) closed 
under the adjoint operation, we define its commutant as 
the von Neumann algebra 

M.' = (x g B(H) VM G M [M, X] = o} . (8) 

Then the Bicommutant Theorem Jl6j states that M." = 
(A4')' is the smallest von Neumann algebra containing 
M.. In particular, when M. is already an algebra, M." is 
the weak closure of M. . Von Neumann algebras are char- 
acterized by the property M" = A4. A von Neumann 
algebra M. with the property that its only elements com- 
muting with all others are the multiples of the identity 
(i.e., M 1 n M" = CI) is called a factor. 

It might seem that the two ways out of the NoGo- 
Theorem indicated at the end of the previous section are 
opposite to each other, but in fact they are closely re- 
lated. For if a; is a state on a C*-algebra C D A U B, we 
can associate with it a Hilbert space H UJl a representa- 
tion 7r w : C — ► B(H), and a unit vector £1 G H^, such that 
uo(C) — (n,-K u (C)ST), and such that the vectors 7r w (C)f2 
are dense in H^. This is called the Gelfand-Naimark- 
Segal (GNS)-construction [jl3| . Clearly, the given state lo 
is given by a density operator (namely |f2) (il|) in this new 
representation and the algebra can naturally be extended 
to the weak closure tt u (C)". The commutativity of two 
subalgebras is preserved by the weak closure, so the nor- 
mal state |f2)(f2|, and the two commuting von Neumann 
subalgebas tt lj (A)" and tt^B)" are again a bipartite sys- 
tem, which describes essentially the same situation. The 
only difference is that some additional idealized observ- 
ables arise from the weak closure operations, and that 
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some observables in C (those with C > but oj(C) = 0) 
are represented by zero in 7r w . 

We remark that von Neumann's incomplete infinite 
tensor product of Hilbert spaces can be seen as a spe- 
cial case of the GNS-construction: The infinite tensor 
product of C*-algebras 0, Ai is well-defined (see (ll| Sec 
2.6] for precise conditions), essentially by taking the norm 
completion of the algebra of local observables ^ i Ai, with 
all but finitely many factors Ai £ Ai equal to lj. On 
this algebra the infinite tensor product of states is well- 
defined, and we get the incomplete tensor product as the 
GNS-Hilbert space of the algebra 0^ BiTLi) with respect 
to a the pure product state defined by the reference vec- 
tors Xi- 



V. VON NEUMANN ALGEBRAS WITH 
MAXIMAL ENTANGLEMENT 

A. Characterization and basic properties 

Let us analyze the example given in the last section: 
the bipartite state obtained from the incomplete tensor 
product of singlets in Hoooo- We take as Alice's observ- 
able algebra A the von Neumann algebra generated by 
all local Alice operators (and analogously for Boh). The 
bipartite state on these algebras, given by the reference 
vector ® ^ Xi > then has the following properties 

ME 1 A and B together generate B(TL) as a von Neu- 
mann algebra, so there are no additional observ- 
ables of the system beyond those measurable by 
Alice and Bob. 

ME 2 A and B are maximal with respect to mutual com- 
mutativity. (i.e., A — B' and B = A') 

ME 3 The overall state is pure, i.e., given by a vector 

ME 4 The restriction of this state to either subsystem is 
a trace, so uj(AiA 2 ) = lu(A 2 Ai), for A\,A 2 € A. 

ME 5 A is hyperfinite, i.e., it is the weak closure of an 
increasing family of finite dimensional algebras. 

These properties, except perhaps ME 2 (see |l7|) are 
immediately clear from the construction, and the prop- 
erties of the respective local observables. They are also 
true for finite dimensional maximally entangled states on 
H = H A <£> H B , A = B{U A ) ® 1, and B = 1 <g> B(H B )- 
This justifies calling this particular bipartite system max- 
imally entangled, as well. 

There are many free parameters in this construction. 
For example, we could take arbitrary dimensions di < oo 
for the i th pair. However, all these possibilities lead to 
the same maximally entangled system: 

Theorem 2 All bipartite states on infinite dimensional 
systems satisfying conditions ME 1 - ME 5 above are 
unitarily isomorphic. 



Proof. (Sketch). We first remark that A has to be a 
factor, i.e., A H A' = CI. Indeed, using ME 1 and ME 
2, we get A n A' = B' n A' = (B U A)' = B(H)' = CI. 

Now consider the support projection S G A of the 
restriction of the state to A. Thus 1 — S is the largest 
projection in A with vanishing expectation. Suppose that 
this projection does not lie in the center of A, i.e., there is 
an A e A such that AS ^ SA. Let X = (t-S)AS, which 
must then be nonzero, as AS-SA = ((1- S) + S)(AS - 
SA) = X - SA(t - S). Then using the trace property 
we get uj{X*X) = lj(XX*) < \\A\\ 2 lu{1 - S) = 0, which 
implies that the support projection of X* X has vanishing 
expectation. But since X*X < \\A\\ 2 S, this contradicts 
the maximality of (1 — S). It follows that S lies in the 
center of A and that S = 1, because A is a factor. To 
summarize this argument, u> must be faithful, in the sense 
that A <E A, A > 0, and u(A) = imply A = 0. 

Now consider the subspace spanned by all vectors of 
the form AO,, with A £ A. This subspace is invariant 
under A, so its orthogonal projection is in A' — B. But 
since (1 — P) obviously has vanishing expectation, the 
previous arguments, applied to B imply that P = 1. This 
is to say that AQ is dense in Ti or, in the jargon of oper- 
ator algebras, that f2 is cyclic for A. Thus Ti. is unitarily 
equivalent to the GNS-Hilbert space of u> restricted to A, 
and the form of B = A 1 is completely determined by this 
statement. Now a factor admits at most one trace state, 
so w is uniquely determined by the isomorphism type of 
A as a von Neumann algebra, and it remains to show 
that A is uniquely determined by the above conditions. 

A is a factor admitting a faithful normal trace state, so 
it is a "type Ili-factor" in von Neumann's classification. 
It is also hyperfinite, so we can invoke a deep result of 
Alain Connes stating that such a factor is uniquely 
determined up to isomorphism. □ 

For the rest of this section we will study further prop- 
erties of this unique maximally entangled state of infi- 
nite entanglement. The items ME 6, ME 7 below are 
clear from the above proof. ME 8 follows by splitting 
the infinite tensor product either into a finite product 
and an infinite tail, or into factors with even and odd la- 
bels, respectively. ME 9 - ME II are treated in separate 
subsections as indicated. 

ME 6 A and B are factors: A(~)A' = CI. 
ME 7 AVt and BVL are dense in H. 

ME 8 The state contains infinite one-shot entanglement, 
which is not diminished by extracting entangle- 
ment. Moreover, it is unitarily isomorphic to two 
copies of itself. 

ME 9 Every density operator on TL maxim ally violates 
the Bell-CHSH inequality (see Section |VB|). 



ME 10 The gener alize d Schmidt spectrum of fl is flat 
(see Section |V C|). 
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ME 11 Every A G A is complete ly co rrelated with a 
"double" B eB. (see Section |vd| ). 



B. Characterization by violations of Bell's 
inequalities 

If we look at systems consisting of two qubits, maxi- 
mally entangled states can be characterized in terms of 
maximal violations of Bell- inequalities. It is natural to 
ask, whether something similar holds for the infinite di- 
mensional setting introduced in Section |v|. To answer 
this question consider again a bipartite state w on an al- 
gebra containing two mutually commuting algebras A, B 
describing Alice's and Bob's observables, respectively. 
We define the Bell correlations with respect to A and 
B in to as 

P(lu) = i supc;(A 1 (B 1 + B 2 ) + A 2 {B 1 - B 2 )), (9) 

where the supremum is taken over all selfadjoint Ai G A, 
Bj G B satisfying -1 < A { < 1, -1 < Bj < 1, for i,j = 
1, 2. In other words A\,A 2 and Bi,B 2 are (appropriately 
bounded) observables measurable by Alice respectively 
Bob. Of course, a classically correlated (separable) state, 
or any other state consistent with a local hidden variable 
model 1 19 satisfies the Bell-CHSH-inequality (3(lu) < 1. 

Exactly as in the standard case, we can show Cirelson's 
inequality (2^, ^l], ^2| bounding the quantum violations 
of the inequality as 



< V2. 



If the upper bound vz is attained we speak of a maximal 
violation of Bell's inequality. 

It is clear that the maximally entangled state described 
above does saturate this bound: In the infinite tensor 
product construction of TL = 7ioooo we only need to take 
observables Ai,Bi from the first tensor factor. But we 
could also have chosen similar observables Ai.k, Bi.k (i — 
1, 2) for the fc th qubit pair. Let us denote by 



T fc = A ltk (B lik + B 2<k ) + A 2 , k {B lt 



B 



2.k 



(11) 



the "test operator" for the k th qubit pair, whose expec- 
tation enters the Bell-CHSH-inequality. Then for a dense 
set of vectors </> G Tt, namely for those differing from the 
reference vector in only finitely many positions, we get 
(</>, T k <f>) = \J~2 for all sufficiently large k. Since the norms 
llXfcH are uniformly bounded, a simple 3e-argument shows 
that lim / ! c _ +oo (0, T k 4>) = \/2 for all 4> G 7ioooo- By taking 
mixtures we find 



lim tr(pT fc ) = V2 



(12) 



for all density operators p on 7i oooo ■ 

This property is clearly impossible in the finite dimen- 
sional case: any product state would violate it. This 
clarifies the statement in Section IV C that TLoooo is in 



no way a tensor product of Hilbert spaces for Alice and 
Bob. Of course, we can simply define a product state 
on the algebra of local operators, and then extend it by 
the Hahn-Banach Theorem to all operators on BCHoooo)- 
However, just as the reference state of infinitely many 
singlets is a singular state on B{TL q®'Hoo), any product 
state will necessarily be singular on B{TL 0000 ). 

It is interesting that bipartite states with property ( |l2] ) 
naturally arise in quantum field theory, with A and B the 
algebras of observables measurable in two causally dis- 
joint (but tangent) spacetime regions. This is true under 
axiomatic assumptions on the structure of local algebras, 
believed to hold in any free or interacting theory. The 
only thing that enters is indeed the structure of the local 
von Neumann algebras, as shown by the following Theo- 
rem || ||, ^l| . Again the maximally entangled state plays 
a key role. 

Theorem 3 (Q) Let A,Bc B(TL) be mutually commut- 
ing von Neumann algebras acting on a separable Hilbert 
space Tl. Then the following are equivalent: 

(i) For some density operator p, which has no zero eigen- 

values, we have (3{p) = 

(ii) For every density operator p on H we have /3(p) = 

V2. 

(iii) There is a set T k of test operators formed from A 

and B such that fiL!\ ) holds for all density operators 
P- 



(10) (i y ) There is a unitary isomorphism under which 



A = Ax® A, 
B = Bi®B, 

Ai,Bi C B(Hoooo) are the algebras of Theorem^, 
and A,B C B{TL) are other von Neumann algebras. 

In other words, the maximal violation of Bell's inequal- 
ities for all normal states implies that the bipartite sys- 
tem is precisely the maximal entangled state, plus some 
additional degrees of freedom (A, B) , which do not con- 
tribute to the violation of Bell inequalities. 

C. Schmidt decomposition and modular theory 

The Schmidt decomposition is a key technique for an- 
alyzing bipartite pure states in the standard framework. 
It represents an arbitrary vector f2 G Ha €5 H-B as 



(13) 



where the c a > are positive constants, and {e a } C Ha 
and {f a } C Hb are orthonormal systems. 
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Its analog in the context of von Neumann algebras is a 
highly developed theory with many applications in quan- 
tum field theory and statistical mechanics, known as the 
modular theory of Tomita and Takesaki [P3| . We recom- 
mend Chapter 2.5 in Jl3| for an excellent exposition, and 
only outline some ideas and indicate the connection to 
the Schmidt decomposition. 

Throughout this subsection, we will assume that 
A, B C B(H) are von Neumann algebras, and ft S TL 
is a unit vector, such that the properties ME 2, ME 
3, and ME 7 of Section VA hold. As in the case of 



the usual Schmidt decomposition the essential informa- 
tion is already contained in the restriction of the given 
state to the subalgebra A, i.e., by the linear functional 
u>(A) = (fi, ACl). Indeed, the Hilbert space and the cyclic 
vector il (cf. ME 7) satisfy precisely the conditions for 
the GNS-representation, which is unique up to unitary 
equivalence. Moreover, condition ME 2 fixes B as the 
commutant algebra. 

However, since A often does not admit a trace, we 
cannot represent w by a density operator, and therefore 
we cannot use the spectrum of the density operator to 
characterize u>. Surprisingly, it is equilibrium statisti- 
cal mechanics, which provides the notion to generalize. 
In the finite dimensional context, we can consider every 
density operator as a canonical equilibrium state, and de- 
termine from it the Hamiltonian of the system. This in 
turn defines a time evolution. Note that the Hamiltonian 
is only defined up to a constant, so we cannot expect to 
reconstruct the eigenvalues of H , but only the spectrum 
of the Liouville operator a i— > i[cr,H], which generates 
the dynamics on density operators, and has eigenvalues 
i(E n — E m ), when the E n are the eigenvalues of H . The 
connection between the time evolutions and equilibrium 
states makes sense also for von Neumann algebras, and 
can be seen as the physical interpretation of modular 
theory @. 

We begin the outline of this theory with the anti-linear 
operator S on TL by 



S(A(l) = A*n, A e A. 



(14) 



It turns out to be closable, and we denote its closure by 
the same letter. As a closed operator S admits a polar 
decomposition 

5=JA 1 / 2 , (15) 

which defines the anti-unitary modular conjugation J and 
the positive modular operator A. 

Let us calculate A in the standard situation, where 
H = JC®]C, and A = BQC)®t respectively B = 1&BQC), 
and f2 is in Schmidt form (|l3|). Due to assumption ME 
7 (cyclicity) , the orthonormal systems e a and f a have to 
be even complete (i.e., bases). Now consider (|l^) with 
A = (|e i o)(e 7 |) ® 1, which becomes 



Sfaep ® / 7 ) = cpe-f ® fp, 



(16) 



from which we readily get 



where p = J2 a c\ |e Q )(e Q | is the reduced density operator, 
F<pi ® 02 = 02 ® <f>i is the flip operator and <d denotes 
complex conjugation in the e„ basis. The time evolution 
with Hamiltonian H = — log/? + cl, for which to is now 
the equilibrium state with unit temperature, is then given 
by £ t (A) <g> 1 = A U (A® t)A- lt . 

In the case of general von Neumann algebras, the spec- 
trum of A need no longer be discrete, and it can be a 
general positive, but unbounded selfadjoint operator. It 
turns out that A** still defines a time evolution on the al- 
gebra A, the so-called modular evolution The equilibrium 
condition cannot be written directly in the Gibbs form 
p oc exp(— H), since there is no density matrix any more, 
but has to be replaced by the so-called KMS-condition, a 
boundary condition for the analytic continuation of cor- 
relation functions |^| which links the modular evo- 
lution to the state. 

In the standard situation, the eigenvalue 1 of A plays 
a special role, because it points to degeneracies in the 
Schmidt spectrum. In the extreme case of a maximally 
entangled state all c Q are equal, and A = 1 or, equiva- 
lently, S is anti-unitary. This characterization of maxi- 
mal entanglement carries over to the von Neumann alge- 
bra case: S is anti- unitary if and only if for all A\, A 2 G A 



{Vt,A x A 2 £l) 



(A* 1 n,A 2 n) 

{A^n,AiQ) 



(SA&sAzn) 
(n,A 2 A 1 n). 



A l/2 



>p~ 1/2 , and J = F(e®G), (17) 



This is precisely the trace property ME 4. 



D. Characterization by the EPR-doubles property 

In the original EPR-argument it is crucial that certain 
observables of Alice and Bob are perfectly correlated, so 
that Alice can find the values of observables on Bob's side 
with certainty, without Bob having to carry out this mea- 
surement. An approach to studying such correlations was 
proposed recently by Arens and Varadarajan J25|. The 
basic idea, stripped of some measure theoretic overhead, 
and extended to the more general bipartite systems con- 
sidered here p6| , rests on the following definition. Let 
A, B be commuting observable algebras and ui a state on 
an algebra containing both A and B. Then we say that 
an element B 6 B is an EPR- double of A € A, or that A 
and B are doubles (of each other) if 

w({A* - B*){A- B)) =oj({A-B)(A*-B*)) =0. (18) 

Of course, when A and B are hermitian, the two expres- 
sions coincide, and in this case there is a simple inter- 
pretation of equation (|l8l). Since A and B commute, we 
can consider their joint distribution (measuring the joint 
spectral resolution of A and B). Then (A — B) 2 is a posi- 
tive quantity, which has vanishing expectation if and only 
if the joint distribution is concentrated on the diagonal, 
i.e., if the measured values coincide with probability one. 

Basic properties are summarized in the following 
Lemma. 
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Lemma 4 Let lo be a state on a C*-algebra containing 
commuting subalgebras A and B. Then 

(i) A and B are doubles iff for all C in the ambient 

observable algebra we have ui(AC) = lo(BC) and 
lo{CA)=lo(CB). 

(ii) IfAi,A2 have doubles B\, B2, then A\, A\+ A2, and 

A1A2 have doubles Bl, Bi + B2, and B2B1, respec- 
tively. 

(iii) When A and B are normal (A A* = A* A), and 

doubles of each other, then so are f(A) and f(B), 
where f is any continuous complex valued function 
on the spectrum of A and B, evaluated in the func- 
tional calculus. 

(iv) When A and B are von Neumann algebras, and lo is 

a normal state, and observables A n with doubles B n 
converge in weak* -topology to A, then every cluster 
point of the sequence B n is a double of A. 

(v) Suppose that lo restricted to B is faithful (i.e., B 3 

B > and lo(B) = imply B = 0). Then every 
A iE A admits at most one double. 

Proof. (i) One direction is obvious by setting C — 
A* — B* . The other direction follows from the Schwartz 
inequality \lo(X*Y)\ 2 < lo(X*X)lo(Y*Y). 

The remaining items follow directly from (i) . (iii) is ob- 
vious from (ii) for polynomials in A and A* , and extends 
to continuous functions by taking norm limits on the 
polynomial approximations to / provided by the Stone- 
Weierstrafi approximation theorem. For (iv) one has to 
use the weak*-continuity of the product in each factor 
separately (see e.g. j27| Theorem 1.7.8]). □ 

In the situation we have assumed for modular theory, 
we can give a detailed characterization of the elements 
admitting a double: 

Proposition 5 Suppose A and B = A' are von Neu- 
mann algebras on a Hilbert space hi, and the state lo is 
given by a vector n E hi, which is cyclic for both A and 
B. Then for every A E A the following conditions are 
equivalent: 

(i) A has an EPR- double B E B. 

(ii) A is in the centralizer of the restricted state, i.e., 

u(AAi) = u(AxA) for all A x E A. 



(iii) A is invariant under the 
A it AA~ it = A for all t E M. 



modular evolution 



In this case the double is given by B = J A* J. 

Proof. (i)=>(ii) When A has a double B, we get 
w{AAi) = to(BAi) = oj(AiB) = uj(AiA) for all A\ in 
the ambient observable algebra. 

(ii)<J=>(iii) This is a standard result (see, e.g., [psl Prop. 
15.1.7]). 



(iii)=>(i) Since A lt n = O, (iii) implies A lt AfL = Aft, 
so Afl is an eigenvector for eigenvalue 1 of the unitary 
A lt and AAfl = Ati. By the same token, AA*fl = A* CI. 
We claim that in that case B = J A* J E B is a double 
of A in B: We have Bfl = J A* JO, = JA*Q = JSAtt = 
A An = An and, similarly, B*n = A*n. From this (i) 
follows immediately. 

The formula for B was established in the last part of 
the proof. Uniqueness follows from Lemma ||. □ 

Two special cases are of interest. On the one hand, in 
the standard case of a pure bipartite state we get a com- 
plete characterization of the observables which posses a 
double: they are exactly the ones commuting with the 
reduced density operator plfl . On the other hand, we 
can ask under what circumstances all A E A admit a 
double. Clearly, this is the case when the centralizer in 
(ii) of the Proposition is all of A, i.e., if and only if the 
restricted state is a trace. Again this characterizes the 
everybody's maximally entangled states on finite dimen- 
sional algebras, and the unique infinite dimensional one 
for hyperfinite von Neumann algebras. 



VI. THE ORIGINAL EPR STATE 

In their famous 1935 paper M Einstein, Podolsky and 
Rosen studied two quantum particles with perfectly cor- 
related momenta and perfectly anticorrelated positions. 
It is immediately clear that such a state does not ex- 
ist in the standard framework of Hilbert space theory: 
the difference of the positions is a self-adjoint operator 
with purely absolutely continuous spectrum, so what- 
ever density matrix we choose, the probability distri- 
bution of this quantity will have a probability density 
with respect to Lebesgue measure, and cannot be con- 
centrated on a single point. Consequently, the wave func- 
tion written in Q is a pretty wild object. Essentially it 
is ^(a^i, X2) — cS(xi — X2 +a), with the Dirac delta func- 
tion, and c a "normalization factor" which must vanish, 
because the normalization integral for the delta function 
is undefined, but infinite if anything. 

How could such a profound physical argument be based 
on such an ill-defined object? The answer is probably 
that the authors were completely aware that they were 
really talking about a limiting situation of more and more 
sharply peaked wave functions. We could model them by 
a sequence of more and more hig hly sq ueezed two mode 
Gaussian states (cf. Subsection VI E[ ), or some other 
sequence representation of the delta function. The key 
point is that the main argument does not depend on the 
particular approximating sequence. But then we should 
also be able to discuss the limiting situation directly in 
a rigorous way, and extract precisely what is common to 
all approximations of the EPR state. 
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A. Definition 

In this section we consider a family of singular states, 
which describes quite well what Einstein Podolsky and 
Rosen may have had in mind. Throughout we assume we 
are in the usual Hilbert space Ti — C 2 (M. 2 ) for describing 
two canonical degrees of freedom, with position and mo- 
mentum operators Q±, Q2, Pi, P%- The basic observation 
is that the operators P\ + Pi and Qi — Q2 commute as 
a consequence of the Heisenberg commutation relations. 
Therefore we can evaluate in the functional calculus (i.e., 
using a joint spectral resolution) any function of the form 
g(Pi + P2, Qi — Q2), where g : M 2 — > C is an arbitrary 
bounded continuous function. We define an EPR-state 
as any state u> such that 

w(ff(Pi+P 2 ,Qi-Q 2 )) =9(0,a), (19) 

where a is the fixed distance between the particles. Sev- 
eral comments are in order. First of all, if we take 
any sequence of vectors to "approximate" the EPR wave 
function (and adjust normalization on the way), weak*- 
clustcr points of the corresponding sequence of pure 
states exist by compactness of the state space, and all 
these will be EPR states in the sense of our definition. 
Secondly, condition (^) does not fix to uniquely. Indeed, 
different approximating sequences may lead to different 
lo. Even for a fixed approximating sequence it is rarely 
the case that the expectation values of all bounded oper- 
ators converge, so the sequence will have many different 
cluster points. Thirdly, the existence of EPR states can 
also be seen more directly: the algebra of bounded con- 
tinuous functions on R 2 is faithfully represented in B(TL) 
(i.e., g(Pi + P2,Qi — Q2) = only when g is the zero 
function). On that algebra the point evaluation at (0, a) 
is a well defined state, so any Hahn-Banach extension of 
this state to all of B(7i) will be an EPR state 4 . 

In our further analysis we will only look at properties 
which are common to all EPR states, and which arc hence 
independent of any choice of approximating sequences. 
The basic technique for extracting such properties from 
(|l9| ) is to use positivity of lo in the form of the Schwartz 
inequality \uj{A* B)\ < lo(A* A)lo(B* B). For example, we 
get 

uo(Xg) = u(gX^=g(0,a)uj(X) , (20) 



4 The reason for defining EPR-states with respect to continuous 
functions of Pi + P2 and Qi — Q2 rather than, say, measurable 
functions, is that we need faithfulness. The functional calculus 
is well defined also for measurable functions, but some functions 
will evaluate to zero. In particular, for the function g(p, x) = 1 
for x = a and p = 0, but g(p, x) = for all other points, we get 
g(Pl + P2,Qi — Q2) = 0, because the joint spectrum of these 
operators is purely absolutely continuous. Hence condition (|is|), 
extended to measurable functions would require the expectation 
of the zero operator to be 1. 



where g is shorthand for g(P\ + P2, Qi — Q2) for some 
bounded continuous function g, and X G B{TL) is an 
arbitrary bounded operator. This is shown by taking 
A = X* and B = (g - g(0,a)t) (or A = (g- g(0,a)l) 
and B = X) in the Schwartz inequality. 

B. Restriction to the CCR-algebra 

Next we consider the expectations of Weyl operators 

W(d,&,»7i,»b) = S^+^-^-mQi) 

= e Ki-P-v-Q) . (21) 

Obviously, if £1 = £2 and 771 = —772, which we will abbre- 
viate as (£, ff) € S, we have W(£, ff) =g for a uniformly 
continuous g, so ( |l9| ) determines the expectation. Com- 
bining it with Equation (^(i| ) we get: 

for {lff)eS. (22) 

In particular, the state is invariant under all phase space 
translations by vectors in S. 

This is already sufficient to conclude that the state is 
purely singular, i.e., that w(K) = for every compact op- 
erator, and in particular for all finite dimensional projec- 
tions. An even stronger statement is that the restrictions 
to Alice's and Bob's subsystem are purely singular. 

Lemma 6 For any EPR state, and any compact operator 
K, lo(K <g> 1) = 0. 

Proof. Indeed the restricted state is invariant under all 
phase space translations, since we can extend W(£, 77) 
to a Weyl operator of the total system, i.e., W'(^,7y) = 
W(£, e, V, -V) = W(f, ry)®W(£, -ry), with (£, £, 77, -rf) e 
S, and 

^((W^AW^t?)*)®!) (23) 

= C(W'(^7 7 )(A®1)W'(^,7 7 )*)) . 

Now consider a unit vector \ with bounded sup- 
port in position space, and let K = |x)(x| be the 
corresponding one-dimensional projection. Then suf- 
ficiently widely space translates W(n£o,0)x are or- 
thogonal, and hence, for all N, the operator Kn = 
T,n=i W« o ,0)iTW« o ,0*) is bounded by 1. Hence 
Nlo(K) = uj(K n ) < oj(t) = 1, and oj(K) = 0. Since vec- 
tors of compact support are norm dense in Hilbert space, 
the conclusion holds for arbitrary □ 

For other Weyl operators we get the expectations from 
the Weyl commutation relations 

WCf, ff)W(i' , rf') = e'f W(f + f , ff+ff') , 

with a = £ • ff' - I*' • 77 . (24) 
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This is just a form of the Heisenberg commutation rela- 
tions. Now S is a so-called maximal isotropic subspace 
of phase space, which is to say that the commutation 
phase a vanishes for (£, ff), (£',ff') £ S, and no subspace 
of phase space strictly including S has the same property. 

For a point (£, ff) in phase space, which does not 
belong to S, we can find some vector (£',f)') £ S 
such that the commutation phase e 1<T ^ 1 is non triv- 
ial. Combining the Weyl relations ( pjj ) with the invari- 
ance © gives w(W(£t7)) = w(W(f , 77 ')W(£ 77)) = 
e iff w(W(£»flW(f / ,Tf')) = e i<T w(W(£ rf)) which implies 
that the expectation values 



w(W(£,»fl)=0 for (£,77)^ 



(25) 



must vanish. With equations ( ]22] , |25| ) we have a complete 
characterization of the state u> restricted to the "CCR- 
algebra" , which is just the C*-algebra generated by the 
Weyl operators. Since this is a well-studied object, one 
might make these equations the starting point of an in- 
vestigation of EPR states. However, one can see that 
( |l9| ) is strictly stronger: there are states which look like 
u> on the CCR-algebra, but which give an expectation in 
( |l9| ) corresponding to a limit of states going to infinity 
instead of going to zero. 



C. EPR-correlations 

How about the correlation property, which is so impor- 
tant in the EPR-argument? The best way to show this is 
the 'double' formalism of Section VD| in which we denote 
by Z the norm closed subalgebra of operators on £ 2 (R) 
generated by all operators of the form f(£P+rjQ), where 
/ : M — > C is an arbitrary uniformly continuous function 
evaluated in the functional calculus on a real linear com- 
bination £P + r\Q of position and momentum 5 . This 
algebra is fairly large: it contains many observables of 
interest, in particular all Weyl operators and all compact 
operators. It is closed under phase space translations, 
and these act continuously in the sense that, for Z £ Z, 
HW^ZW^r;)* _ Z \\ -+ as (£,77) -» 6 . 

Theorem 7 All operators of the form Z®\ with Z £ Z 
have doubles in the sense of equation (18). Moreover, the 
double of Z®\ is \®Z T , where Z T denotes the transpose 
( adjoint followed by complex conjugation ) in the position 
representation. 



Proof. We only have to show that for f(£P + rjQ) <g> 1 
f($Pi + vQi) we get the double 1 ® /(-£P + rjQ) 



5 The same type of operators, although motivated by a different 
argument already appears in [ft 

lis continuity is crucial in the correspondence theory set out in 
We where not able to prove the analogue of Theorem M by 
only assuming this continuity. 



/(— ^2+77*32), when /, £, and 77 are as in the definition of 
Z. By the general properties of the double construction 
this will then automatically extend to operator products 
and norm limits. 

Fix e > 0. Since / is uniformly continuous, there is 
some S > such that \ f(x) — f(y)\ < e whenever \x — y\ < 
5. Now pick a continuous function h : K — > [0, 1] C K 
such that h(0) = 1, h(t) = for \t\ > 5. We consider the 
operator 



M 



(/(CA+77Q1) 

F^Pi + nQu- 



-/(- 
Pi)- 
ip 2 



-eP 2 +77Q 2 )) x 

■v(Qi-Q2)) 

-77Q2) 



where F(x,y) — (/(x) — f(y))h(x — y) and this function 
is evaluated in the functional calculus of the commuting 
selfadjoint operators (£Pi +T)Qx) and (— £P2+r)Q2) ■ But 
the real valued function F satisfies \F(x,y)\ < e for all 
(x,y): when \x — y\ > 5 the /i-factor vanishes, and on the 
strip \x — y\ < 8 we have \f(x) — f(y)\ < s. Therefore 
||M|| < e. Let X be an arbitrary operator. Then 

H[f(tPi + VQi) - !(-iP2 + rtQi)]X) I 
= \lu(MX)\ < \\M\\ \\X\\ < e\\X\\ . 

Here we have added a factor h(£(Pi+P2)+ri(Qi — Q2)) at 
the second equality sign, which we may because of (20), 
and because h is a function of the appropriate operators, 
which is = 1 at the origin. Since this estimate holds for 
any e, we conclude that the first relation in Lemma [I|.l 
holds. The argument for the second relation is completely 
analogous. □ 



D. Infinite one-shot entanglement 

In order to show that the EPR state is indeed highly 
entangled, let us verify that it contains infinite one-shot 
entanglement in the sense forbidden by Theorem [j]. The 
local operations needed to extract a <i-dimensional sys- 
tem will be simply the restriction to a subalgebra. In 
other words, we will construct subalgebras Ad C A and 
Bd C B such that the state u> restricted to Ad <£> Bd will 
be a maximally entangled pure state of ci-dimensional 
systems. 

The matrix algebras Ad,Bd are best seen to be gen- 
erated by Weyl operators, satisfying a discrete version 
of the canonical commutation relations (|2~4|), with the 
addition operation on the right hand side replaced by 
the addition in a finite group. Let Zd denote the cyclic 
group of integers modulo d. With the canonical basis 
\k, I), k,£ £ Id w e introduce the Weyl operators 

w(ni,mi,n 2 ,m 2 )|fc,^) (26) 

= £Mk-m 1 )+n 2 (i-m 2 )\ k _ mi j_ rn2 ^ 

where £ = exp(2iri/d) is the dth root of unity. These 
are a basis of the vector space B(C d ® C d ), which 
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shows that this algebra is generated by the four unitaries 
u x = w(l, 0, 0, 0), v l = w(0, 1, 0, 0), u 2 = w(0, 0, 1, 0) and 
v 2 = w(0, 0, 0, 1). They are defined algebraically by the 
relations v k Uk — C u k v k ,k — 1,2, and uf = u\ — vf — 
i>2 = 1. The one dimensional projection onto the stan- 
dard maximally entangled vector fl = dr 1 ^ 2 \kk) can 
be expressed in the basis (|2^) as 

= ^E^Tc^r, (27) 



which will be useful for computing fidelity. 

In order to define the subalgebras extracting the de- 
sired entanglement we first define operators Ui , V\ in Al- 
ice's subalgebra and U 2 ,V 2 m Bob's, which satisfy the 
above relations and hence generate two copies of the d x d 
matrices. It is easy to satisfy the commutation relations 
VkUk = (UkVk, by taking appropriate Weyl operators, 
say 



Ui 



u 2 



and V k = e* Pk (28) 



with £ = 2n/d. The tilde indicates that these are 
not quite the operators yet we are looking for, because 
they do not satisfy the periodicity relations: Uf — 
exp(idQi) 1, and similarly for U 2 and Vk- We will 
denote by A the C*-algebra, generated by the operators 
U\,V\ (p8|). The algebra B is constructed analogously. 
Then by virtue of the commutation relations Uf and Vf 
commute with all other elements of A, i.e., they belong to 
the center Ca C A, which represents the classical vari- 
ables of the system. In the same manner, U% and V 2 
generate the center Cb of Bob's algebra B 7 . 

If we take any continuous function (in the functional 
calculus) of a hermitian or unitary element of Ca, it will 
still be in Ca- If we take a measurable (possibly discon- 
tinuous) function the result may fail to be in Ca, but it 
still commutes with all elements of A (and analogously 
for Bob's algebras). In particular, we construct the op- 
erators 



u k = (u d k ) 1/d , 



(29) 



where the dth root of numbers on the unit circle is taken 
with a branch cut on the negative real axis. This branch 
cut makes the function discontinuous, and also makes 
this odd-looking combination very different from Uk ■ We 
now define V k analogously, and set 



U k = U7 x U k and V k = Vr x V k 



(30) 



for k = 1, 2. Then since U k , V k commute with A®B, the 
commutation relations VkUk = QUkVk still hold, but in 



addition we have U% 



1, because U?. 



U d k - 



It remains 



to show that on the finite dimensional algebras gener- 
ated by these operators, the given state is a maximally 
entangled pure state. We will verify this by computing 
the fidelity, i.e., the expectation of the projection (p7H: 



^^(Uiu^nwT 1 ) = 1. 



(31) 



Proof of this equation. We have shown in Section VI C 



that Ui and U 2 are EPR-doubles. This property trans- 
fers to arbitrary continuous functions of U% and U 2 by 
Lemma ^ and uniform approximation of continuous func- 
tions by polynomials. However, because the state ui is not 
normal, it does not transfer automatically to the measur- 
able functional calculus and hence not automatically to 
U\ and U 2 . We claim that this is true nonetheless. 

Denote by ra(z) — z 1 ^ the dth root function with 
the branch cut as described, and let f e be a continuous 
function from the unit circle to the unit interval [0, 1] 
such that f e (z) — 1 except for z in an e-neighborhood of 
2 = —1 in arclength, and such that / e (— 1) = 0. Then the 
function z 1— > f e (z)rd(z) is continuous. Then, since Uf 
and U 2 are doubles, so are f t (U d ), f e (Uf)Ui and their 
counterparts. Note that both of these commute with 
all other operators involved. Hence (using the notation 
\X\ 2 = X*X or \X\ 2 = XX*, which coincide in this case) 

= u>^f e {U*% - MU d )U2\) = 0, (32) 

where the first equality holds by expanding the modulus 
square, and applying the double property of fe(Uf ) where 
appropriate. On the other hand, we have 



7 The C*-algebra A is isomorphic to the continuous sections in an 
C*-algebra bundle over the torus, where each fiber is a copy of 
the algebra Ad- Such a bundle is called trivial, if it is isomorphic 
to the tensor product Ad <8> Ca ■ This would directly give us the 



desired subalgebra J- 
bundle is not trivial 



subalgebra of A- However, this is 
2SL|30|. In order to "trivialize" the bundle, 
we are therefore forced to go beyond norm continuous operations, 
which respect the continuity of bundle sections. Instead we have 
to go to the measurable functional calculus, and introduce an 
operation on the fibers, which depends discontionuously on the 
base point, through the introduction of a branch cut. 



«, [h.-uutf)\u 1 -u 2 



<4u(t-f e (U?) 2 ) <4- 



(33) 



because ||Z7i-Z7 2 || < 2, and < f e {Uf) < 1. For the esti- 
mate we used that fe(z) 2 for all z on the unit circle except 
a section of relative size 2e/(27r), and that the probability 
distribution for the spectrum of Uf is uniform, because 
the expectation of all powers (Uf) = exp(indQ\) van- 
ishes. 
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Adding @ and we find that w(|[/i-£/ 2 | 2 ) < 4e/7r 
for every e, and hence that U\ and J72 are EPR doubles as 
claimed. The proof that V\ and V 2 * are likewise doubles 
(just as Vi and V 2 *) is entirely analogous. Hence U\ and 
U2 as well as Vi and V 2 are also doubles. Applying this 
property in the fidelity expression (|3l|) we find that every 
term has expectation one, so that with the prefactor d~ 2 
the d 2 terms add up to one as claimed. □ 



E. EPR states based on two mode Gaussians 

In this section we will deviate from the announcement 
that we intended to study only such properties of EPR 
states which follow from the definition alone, and are 
hence common to all EPR states. The reason is that there 
is one particular family, which has a lot of additional 
symmetry, and hence more operators admitting doubles, 
than general EPR states. Moreover, it is very well known. 
In fact, most people working in quantum optics probably 
have a very concrete picture of the EPR state, or rather 
of an approximation to this state: since Gaussian states 
play a prominent role in the description of lasers, it is 
natural to consider a Gaussian wave function of the form 



1 / 1 - A 

A exp l-iOTA) 

l + A 
4(1 -A) 

OO 



(?i - 92) 2 

(qi + ?2) 2 



(34) 
(35) 



where e„ denotes the eigenbasis of the harmonic oscil- 
lators Hi = (P 2 + Qi)/2 (i = 1,2). This state is also 
known as the NOPA state, and the parameter A S [0, 1) 
is related to the so-called squeezing parameter r by 
A = tanh(r). Values around r = 5 are considered a 
good experimental achievement |3l} ]. Of course, we are 
interested in the limit r — ► 00, or A — ► 1. 

The A-dependence of the wave function can also be 
written as 

^x(xi,x 2 ) = ^ (xichr] + X2shri, -x 1 sli77 + x 2 ch?7) , (36) 

where the hyperbolic angle 77 is r/2. It is easy to see that 
for any wave function 9 the probability distributions of 
both Qi — Q2 and Pi + P2 scale to a point measures at 
zero. Hence any cluster point of the associated sequence 
of states uj\(X) = (^\, A*f?\) is an EPR state in the sense 
of our definition (with shift parameter a = 0). Note, 
however, that the family itself does not converge to any 
state: it is easy to construct observables X for which the 
expectation oj\(X) remains oscillating between and 1 
as A — > 1. Here, as in the general case, a single state can 
only be obtained by going to a finest subsequence (or by 
taking the limit along an ultrafilter). 



The virtue of the particular family (|35|) is that it has 
especially high symmetry: it is immediately clear that 



{(J{Hi)-f{H 2 ) 



*A =0 



(37) 



for all A, and for all bounded functions / : N — > C of 
the oscillator Hamiltonians Hi,H 2 . This implies that 
f(Hi) and f(H 2 ) are doubles with respect to the state 
lu\ for each A. Clearly, this property remains valid in the 
limit along any subsequence, so all EPR-statcs obtained 
as cluster points of the sequence uj\ also have f(Hi) 
in their algebra of doubles. Consequently, the unitaries 
Uk(t) — exp(itHk) are also doubles of each other, and 
the limiting states are invariant under the time evolution 
Ui2(t) = U\(t) ® U 2 (—t). This is certainly suggestive, 
because oscillator time evolutions have an interpreta- 
tion as linear symplectic transformations on phase space: 
Qk >— > Qk cost ± Pk sini and P% 1— ► TQk sini + Pt cost, 
where the upper sign holds for k = 1 and the lower for 
k = 2. The subspace S from Section VIC is invariant 
under such rotations, and one readily verifies that the 
time evolution U\2{t) takes EPR states into EPR states. 
This certainly implies that by averaging we can generate 
EPR states invariant under this evolution, and we have 
clearly just constructed a family with this invariance. 

As A — * 1, the Schmidt spectrum in ( ^B|) becomes 
"flatter" , which suggests that exchanging some labels n 
should also define a unitary with double. Let p : N — > N 
denote an injective (i.e., one-to-one but not necessarily 
onto map). Then we define an isometry V p by 



with adjoint 



r 1 e p -i (n) if n €p(N) 

" " ~ S if n £ p(N) 



(38) 



(39) 



Let us assume that p has finite distance, i.e., there is a 
constant £ such that \p(n) — n\ < £ for all n € N. We 
claim that in this case V r , 



)t and 1®V* are doubles in all 



EPR states constructed from the sequence (p5|). We show 
this by verifying that the condition holds approximately 
already for finite A. Consider the vector 



(40) 



\/l-A 2 - A p(n) ) e p(n) ® e„ 



n=0 



where in the second summand we changed the summation 
index from n to p(n), automatically omitting all terms 
annihilated by V* according to ( |39"| ) . Since this is a sum 
of orthogonal vectors, we can readily estimate the norm 
by writing (A™ - X p ^) = A"(l - \PW- n ): 

\\A X \\ 2 <max|l-A p(n) -™| 2 < |1 — A ^ | 2 , (41) 

n 

which goes to zero as A — ► 1. Therefore 

cj a (x(F p ®1-1®F;)) = XA\) — > (42) 
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as A — > 1. Hence V p ® 1 and 1 £g> V^,* are doubles in any 
state defined by a limit of u)\ along a subsequence, as 
claimed. 

V p is an isometry but not necessarily unitary. But it 
is effectively unitary under an EPR state: Since V p is 
in the centralizer, we must have — V P V*) ® T) — 

— V*V P ) ® l) = 0, although this operator is non- 
zero. This is in keeping with the general properties of 
EPR states, whose restrictions must be purely singular. 
In fact, (1—VpVp*) is the projection onto those eigenstates 
e„ for which n £ p(N), and this set is finite: it has at 
most t elements 8 . 

It is interesting to note what happens if one tries to re- 
lax the finite distance condition. An extreme case would 
be the two isometries Kven e„ = e2„ and V && e„ — 
e2n+i- These cannot have doubles in any state, because 
the restriction loa of the state to the first factor would 
then have to satisfy 1 = ^(KjvenKren + ^oddKfdd) = 

^(KfvenKven + K*dd^,dd) = 0J A (1 + 1) = 2. On the 

other hand, the norm of Ax no longer goes to zero, and 
we get ||Aa|| 2 — > 1/6 instead. 

To get infinite one-shot entanglement is easier than in 
the case of general EPR states: we can simply combine 
d periodic multiplication operators with d-periodic per- 
mutation operators to construct a finite Weyl-system of 
doubles 9 . In fact there is a very quick way to get high 
fidelity entangled pure states even for A < 1 (sec Q for 
an application to Bell inequality violations). Consider 
the unitary operator Ud ■ H — > H <g> C d given by 



U d e dk+r = (e k (g) el d) ) 



(43) 



for k = 0, 1, . . . and r = 0, 1, . . . , d - 1. Then 

(U d ® U d )^x = ® (44) 

with a A-dependent normalized vector <E C d <g> C d 
proportional to 



M/focVA 



(45) 



Note that the infinite dimensional factor on the right 
hand side of (fl4|) is again a state of the form (^), how- 
ever, a less entangled one with parameter A' = X d < A. 
The second factor, i.e., (|44| ) becomes maximally entan- 
gled in the limit A — > 1. Therefore the unitary (U d <X> Ud) 
splits both Alice's and Bob's subsystem, so that the 
total system is split exactly into a less entangled ver- 
sion of itself and a pure, nearly maximally entangled 
d-dimensional pair. The local operation extracting en- 
tanglement from this state is to discard the infinite di- 
mensional parts. Seen in one of the limit states of the 



8 For any N > t, consider the set {1, . . . , TV}. This has to contain 
at least the images of {1, . . . , N— £}, hence it can contain at most 
I elements not in p(N). 

9 This is probably what the authors of |32| are trying to say. 



family u>\ this is maximally entangled, so equation (g) 
is satisfied with e — 0. Moreover, since the remaining 
system is of exactly the same type, the process can be 
repeated arbitrarily often. 



Counterintuitive properties of the restricted 
states 



Basically, subsection [VI Q shows that the EPR states 
constructed here do satisfy the requirements of the EPR 
argument. However, Einstein, Podolsky and Rosen do 
not consider the measurement of suitable periodic func- 
tions of Qk or Pk but measurements of these quantities 
themselves ||: What do EPR states have to say about 
these? 

Unfortunately, the "values of momentum" found by 
Alice or Bob are not quite what we usually mean by 
"values": they are infinite with probability 1. To see 
this, recall the remark after eq. ( |22| ) that EPR states are 
invariant with respect to phase space translations with 
W(f, rj) with (£,77) e S. Hence 

w(W(Ci, 0, m,0)(A ® 1)W(&, 0, 0)*) 

= to (w(a , a, 771, ® 1) w(a ,z um ,- m y) 

= u(A®l). (46) 

That is, the reduced state is invariant under all phase 
space translations. Now suppose that for some con- 
tinuous function / with compact support we have 
w(/(Qi)) = e 7^ 0. Then we could add many (say 
N) sufficiently widely spaced translates of / to get an 
operator F = £)f f(Qi + x,t) with ||F|| < ||/|| and 
|JVe| = \u{F)\ < 11/11, which implies e = 0. Hence 
for every function with compact support we must have 
w(f(Qi)) = 0. Note that this is possible only for sin- 
gular states, since we can easily construct a sequence of 
compactly supported function increasing to the identity, 
whose u) expectations are all zero, hence fail to converge 
to 1. 

In spite of being infinite, the "measured values" of Al- 
ice and Bob are perfectly correlated, which means that 
we have to distinguish different kinds if infinity. Such 
"kinds of infinity" are the subject of the topological the- 
ory of compactifications ]l5| , |34| . The basic idea is very 
simple: consider some C*^algebra of bounded functions 
on the real line. Then the evaluations of the functions at 
a point, i.e., the functionals x 1— * f[x), are pure states on 
such an algebra, but 'compactness of the state space to- 
gether with the Krem-Milman Theorem (35) dictates that 
there are many more pure states. These additional pure 
states are interpreted as the points at infinity associated 
with the given observable algebra. The set of all pure 
states is called the Gel'fand spectrum of the commutative 
C*-algebrafl5|, Sec. 2. 3. 5], and the algebra is known to be 
isomorphic to the algebra of continuous functions on this 
compact space. For the algebra of all bounded function 
the additional pure states are called free ultrafilters, for 
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the algebra of all continuous bounded functions we get 
the points of the Stone-Cech-compactification, and for 
the algebra of uniformly continuous functions we get a 
still coa rser notion of points at infinity. According to 
Section VI C these are the measured values, which will 



be perfectly correlated between Alice's and Bob's posi- 
tions or momenta. It is not possible to exhibit any such 
value, because proving their mere existence already re- 
quires an argument based on the Axiom of Choice. 

So do we have to be content with the statement that 
the measured values lie "out there on the in finite r anges, 
where the free ultrafilters roam?" Section VI D| shows 



that for many concrete problems, involving not too large 
observable algebras, we can use the perfect correlation 
property quite well. A smaller algebra of observables 
means that many points of Gel'fand spectrum become 
identified, and some of these coarser points may have 
a direct physical interpretation. So the moral is not so 
much that compactification points at infinity are wild, 
pathological objects, but that they describe the way a 
sequence can go to infinity in the finest possible detail, 
which is just much finer that we usually want to know. 
The EPR correlation property holds even for such wild 
"measured values" . 
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